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The masses of the positive parity [70, + ] and [70, 2 + ] non-strange baryons are calcu- 
lated in large N c QCD by considering the most dominant operators in an 1/N C expansion. 
The approach is based on the introduction of an excited core, obtained after the last 
particle (an excited quark) has been removed. Configuration mixing is neglected, for sim- 
plicity. Although being a sub- leading l/N c order, we find that the spin-spin interaction 
plays a dominant role in describing the data. The role of N® operators is also pointed 
' out. We show how the contribution of the linear term in N c , of the spin-spin and of the 
spin-orbit terms vary with the excitation energy. 

O 

c-| ■ 1. Introduction 

The large N c QCD approach suggested by 't Hooft [ [T] and considered in detail by 
Witten [ |2] has become a powerful tool in baryon spectroscopy. The method is based on 
the result that, for Nf flavors, the ground state baryons display an exact SU(2iVy) spin- 
flavor symmetry in the large N c limit of QCD [|3J. It has been applied with great success 
^ to the ground state baryons (N = band), described by the symmetric representation 56 
of SU(6), where Nf = 3 [ E2 EJ El |B1 E| • For excited baryons this symmetry is broken. 
However, based on the observation that excited states can be approximately classified 
as SU(2Nf) multiplets, considerable efforts have previously been made to analyze the 
excited states belonging to the [70, l - ] multiplet (N — 1 band) in the large N c limit [ 

noi nn n2i usi us nni nni nzi [thi nni 1201 , with obvious success. 

The N = 2 band contains the multiplets [56', 0+], [56, 2+], [70, 0+], [70, 2+] and [20, 1+]. 
Among them, the baryons supposed to belong to [56, 2 + ] or [56', + ] have been considered 
so far in Refs. [ HHj and [ |^] . The method of Ref . [ EJ] has recently been extended to 
higher excitations belonging to the [56, 4 + ] multiplet (N = 4 band) [12*2*]. For simplicity, 
in these studies configuration mixing has been neglected. 

The symmetric representation 56 requires a much simpler treatment than the mixed 
representation 70. The main reason is that for the symmetric representation it is not 
necessary to distinguish between excited and core quarks. This is possible because of 
the structure of the wave function. Let us consider the case [56, 2 + ]. For large N c this 
multiplet becomes [[N c ],2 + ] 3 . Its intrinsic wave function (center of mass coordinate 
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removed), written in a harmonic oscillator single particle basis, takes the form 



|[N C ],2 + ) = \l^^\[N c ](0s) N - l (0d)) + ^j^\[N c }(0sf^(0p) 2 }, (1) 



c 



which shows that for N c — > oo one can neglect the second term. This cumbersome term, 
involves an excited core after the removal of the last particle, while in the first term 
the entire excitation is carried by one quark, excited to the 1 = 1 shell. The remaining 
symmetric core of N c — 1 unexcited quarks is simpler than an excited core. Now, if we 
consider the spin-orbit 1-body interaction 

H so = w{r)£-s, (2) 
it is easy to show that 

to I /VIM - / °( N c) if X is mixed symmetric (MS) , . 

K 1 1 ; \ 0(N~ 1 ) if X is symmetric (S) ' K ) 

as proved in Ref. [ES] and, in a different manner, in the Appendix below. In this equation 
\I/ and x are the total and the spin-flavor wave functions of the excited baryon respectively. 
(Note that x an d the spatial wave function of the Appendix have the same permutation 
symmetry, thus the same partition.) With this result, the spin-orbit operator takes the 
simple form used in Ref. [ HH] . 

The situation with the multiplets [70, + ] and [70, 2 + ] is different. This can be il- 
lustrated by writing the orbital part of the intrinsic wave functions following the same 
prescription as above. For £ = one obtains 





N c - 1, 1], 0+> p , A = J-\[N e - 1, l] p , x (0s) N °-\ls)) + J-\[N C - 1, lU(0s) N ^ 2 (0p) 2 ). (4) 



In the first term Is is the first single particle radially excited state with n = 1, i = 
(N = 2n + £). In the second term the two quarks are excited to the p-shell to get N = 2. 
They are coupled to £ = 0. The lower indices p and A are constituent quark model 
notations. They distinguish between states which are antisymmetric and symmetric under 
permutation of the first two particles, respectively. By analogy, for £ = 2 one has 



|[N C - 1, 1], 2+> p , A = y||[iV c - 1, l} p , x (0s) N <-\0d)) + y||[iV c - 1, l] PiX (0s) N ^(0p) 2 ), (5) 



where the two quarks in the p-shell are coupled to £ = 2. We have obtained the expressions 
(jlj) and (jHJ) by using the procedure developed in Ref. [Elj based on generalized Jacobi 
coordinates 4 . The N c = 3 case is demonstrated in Ref. [125]. 

One can see that the coefficients of the linear combinations and (jSJ) are independent 
of iV c which means that both terms have to be considered in the large N c limit. In Eqs. 
(HJ) and (J3J) the first term can be treated as in the [70, 1~] sector, i.e. as an excited quark 

[N c ] is d[Ay = jvT^iVj — i)l • With Nf = 3 and N c = 3 one recovers 56. 

4 For consistency with Ref. we use the harmonic oscillator notation \n£) for single particle states 

everywhere in the text. 
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coupled to a ground state core [ HOI HH [T21 H3 US HH1 UHl HZl HH1 CHI 1201 The second 
term will be treated here as an excited quark coupled to an excited core 5 . 

To our knowledge this is the first attempt to incorporate core excitations in the system 
and treat the mass operator accordingly. As there are similarities with the [70, 1~] multi- 
plet, we view this study as an extension of Refs. [ 1101 [TT1 [T^l 1131 [IS 113 IIHl [IZl EH! UHl 120] . 
In practice we shall combine the techniques of Refs. [E] and [I16j. 

Before ending this section we should recall that in reality the excited states are reso- 
nances, so they have a finite width. This is not taken into account in the present approach, 
similarly to constituent quark models where resonances are treated as bound states. How- 
ever, it is important to notice, that the bound state picture turned out to describe the 
baryon phenomenology satisfactorily, thus being a rather realistic picture of baryon reson- 
nances [120]- The decay widths have been considered, for example, in Refs. [ 1201 [20] - The 
conclusion was that a general large N c analysis does not predict narrow widths, which 
would vanish in the large N c limit. Contrary, generic large N c counting gives widths of 
order N®. According to Ref. [120] the narrowness of the excited states is an artifact of 
simple quarks model assumptions used in the calculations, as it is the case here. 



2. The wave function 



We shall separately discuss the two distinct configurations entering the wave functions 
flU) and ©. 

2.1. The configurations (Os) jVc_1 (ls) and (0s) 7Vc " 1 (0ci) 

In Eqs. (0J) and (0) the first term, containing the configurations (0s) 7Vc_1 (ls) and 
(0s) Nc ~ 1 (0d) respectively, is similar in structure with the wave functions of the [70, 1~] 
multiplet (N c = 3), the difference being that the quark is now excited to the N = 2n+£ = 2 
band (£ = or £ = 2) instead of N = £ = 1. 

The configurations (0s) iVc_1 (ls) and (0s) Nc ~ 1 (0d) are then described by a wave function 
written in the notation of Ref. [ EI] as 




i 

2 

m 2 




I + |;mi,ai) <g> ||;m 2 ,a 2 ) ® \£ q ,m tq ), 



(6) 



where S and I are the total spin and isospin of the baryon, £ = £ q is the angular momentum 
of the excited quark, S c and I c are the core spin and isospin respectively, p = S — I = ±1,0 
and T]/2 = I c — I = ±1/2. The coefficients c p>v will be presented below. The coupling 
coefficients are SU(2) coefficients, as we consider here only non-strange baryons (Nf = 2). 

2.2. The configurations (0s) Nc - 2 (0p) 2 

In order to better understand the 1/N C counting it is useful to rewrite the orbital part of 
the wave function using the fractional parentage technique [EZj by which the last particle 



5 One can also consider two excited quarks and leave the core in the ground state but this situation is 
more complicated. 
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is decoupled. Here we consider the only case treated in the literature, where the mixed 
symmetric state belonging to the representation [iV c — 1,1] has the last particle in the 
second row of the corresponding Young diagram. Then, the second term of (J1J) or (0) 
takes the form 



\[N c -l,l] p>x (0s) N °- 2 (0p) 2 ,l + ) 




%„-i] ((Osfc- 2 (Op)) ^(Op) 
(0s) N ^ 3 (Op) 2 ) <P m (0s), 



(7) 



the decoupling being valid both for £ = and 2. Here all states are normalized. The first 
factor in each term in the right-hand side is a symmetric (N c — l)-particle wave function 
and is a one particle wave function associated to the iV c -th particle. One can see 
that for large N c the coefficient of the first term is 0(1) and of the second 0(N^ 1 ^ 2 ). 
Then, in the large N c limit, one can neglect the second term and take into account only 
the first term in the wave function, where the iV c -tli particle has an £ = 1 excitation. 
A similar decomposition can be made for a symmetric state [N c ] or a mixed symmetric 
state [N c — 1,1] containing one excited quark only (see Appendix). In that case one can 
immediately recover the considerations developed in Ref. [EH] for the spin-orbit or other 
operators. 

When both the decoupled quark and the core are excited, the wave function takes the 
form 



\JJ 3 ;Ih;(£;S = I + p)) = 

E 

m q ,m c ,m£ ,mg ,mi ,7712 ,ai ,Q2 



S 

m s 




( s c 


1 
2 


S ) 


(i. 


1 
2 




\ mi 


m 2 


m s J 




a 2 





E 

n=±l 



I + -; 7711,0!].) ® 



m 2 ,a 2 ) ® \£ q ,m q ) ® \£ c ,m c ), 



(8) 



which is a generalization of (jHJ) still for non-strange baryons. It contains an extra SU(2) 
Clebsch-Gordan coefficient which couples the angular momentum £ q of the excited quark 
to the angular momentum £ c of the excited core. The coefficients c PjV are the same in 
(0) and (JBJ). They are defined in Refs. [HH] and [CHI- We recall that for the symmetric 
(SYM) [N c ] and the mixed (MS) [N c — 1, 1] representations of the permutation group Sn c , 
the non-vanishing coefficients c Pi7? are given by 

(9) 



JvIS 
"0+ 



^SYM 
"0+ 



S[N C + 2(S+1)} 
\ N C (2S + 1) ' 



MS 



(S+l)[N c -2S\ 
\ N C (2S + 1) ' 



^SYM 



-0+ • 



(10) 
(11) 



In Ref. [ EI] they were defined as elements of an orthogonal basis rotation. In terms 
of group theory language these coefficients can be identified with isoscalar factors of the 
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permutation group [ 123 I2E] • They are factors of the Clebsch-Gordan coefficients needed 
in the flavor-space inner product and are related to the position of the last particle in the 
corresponding Young tableaux. 

3. The Mass Operator 

For the [70, £ + ] sector the building blocks which form the operators entering in the mass 
formula consist of the excited core operators t c , S l c , T c a and G* a and the excited quark 
operators t q , s l , t a and g ia . We also introduce the tensor operator 6 

tT = \{t\A)-\k,- 3 i a -h, (12) 

with a = c, b = q or vice versa or a = b = c or a = b = q. For simplicity when a = b, we 
shall use a single index c, for the core, and q for the excited quark so that the operators 
are and respectively. The latter case represents the tensor operator used in the 
analysis of the [70, 1~] multiplet (see e.g. Ref. [ITS]). 



Table 1 

List of operators and the coefficients resulting from the fit with Xdof — 0.83. 



Operator Fitted coef. (MeV) 

Ox = N c l ci = 555 ± 11 

2 = c 2 = 47 ± 100 

3 = JLegWgi*Gi a c 3 = -191 ± 132 

4 = ±(S l c Si + s l S l c ) c 4 = 261 ± 47 



Table 2 

Matrix elements of N. 



Ox 2 3 4 



4 iV[70,2+ 


7 + 
2 


N c 


2 




-6^(^ + 1) 


5 


2 iV[70,2+ 


5 + 
2 


N c 




3) 





ife(iV c + 3) 


4 iV[70, 2+ 


5 + 
2 


N c 


i 

9 






5 
2N C 


4 iV[70,0+ 


3 + 
2 


N c 










5 
2N C 


2 A^[70,2+ 


3 + 
2 


N c 




-3) 





ife(^ + 3) 


4 A^[70, 2+ 


3 + 
2 


N c 


2 
3 







5 
2N C 


2 iV[70,0+ 


1 + 
2 


N c 










ife(iV c + 3) 


4 A^[70,2+ 


1 + 
2 


N c 


-1 




-T2fe(^+ 1 ) 


5 
2N C 



6 The irreducible spherical tensors are denned according to Ref. 
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Table 3 

Matrix elements of A. 





o x 


o 2 


o 3 


o 4 


2 A[70,2+]| + 
2 A[70, 2+]§ + 
2 A [70, 0+]| + 


N c 
N c 
N c 


2 
9 

1 

3 









1 

N c 

1 

N c 

1 

Nr. 



We apply the l/N c counting rules presented in Refs. [ H3] and [ HBJ and use their 
conclusions in selecting the most dominant operators in the practical analysis. For non- 
strange baryons, Table I of Ref. [EI] gives a list of 18 linearly independent operators. If 
the core is excited the number of operators appearing in the mass formula is much larger. 
However due to lack of data, here we have to consider a restricted list. The selection is 
suggested by the conclusion of Ref. [ 13 j, (Nf = 2) and of Ref. [ [TB] (Nf = 3), that 
only a few operators, of some specific structure, bring a dominant contribution to the 
mass. Following the notation of Ref. [ ITB] these are Oi, O2, O3 and O4 exhibited here in 
Table I. The first is the trivial operator of order 0(N C ). The second is the 1-body part 
of the spin-orbit operator of order 0(1) which acts on the excited quark. The third is a 
composite 2-body operator formally of order 0(1) as well. It involves the tensor operator 
(I12j) acting on the excited quark and the SU(6) generators g ia acting on the excited quark 
and G{ a acting on the core. The latter is a coherent operator which introduces an extra 
power N c so that the order of O3 is 0(1), as it can be seen from Table El In order to take 
into account its contribution we have applied the rescaling introduced in Ref. [E] which 
consists in introducing a multiplicative factor of 3. Without this factor the coefficient C3 
becomes too large, as noticed in Ref. [HH] 1 ■ The dynamics of the operator 3 is less 
understood. Previous studies [[T3J[Tn| speculate about its connection to a flavor exchange 
mechanism [ EJUJ EH] related to long distance meson exchange interactions. Finally, the 
last operator is the spin-spin interaction, the only one of order 0(l/iV c ) which we consider 
here. Higher order operators are neglected. Accordingly the mass operator of the [70, l + \ 
multiplet is approximated by 

4 

M [7CM+] =E C ^ > ( 13 ) 
1=1 

where the coefficients q have to be found in a numerical fit to the available data, as 
described below. The diagonal matrix elements of the operators Oi (i = 1, 4) are given 
in Tables El and El They have been obtained from the wave functions (0) and ©. Each 
matrix element contains the additional contribution of the two terms of the wave function. 
We derived general analytic formulae for each case in terms of N c . For the first term which 
has a ground state core, we used Eq. (0) to recover the matrix elements of Appendix A 
of Ref. [HB]- For the second term, using Eq. (|SJ), we have derived analytic expressions 
which generalize those of Ref. [ ^] to an excited core with angular momentum l c 7^ 0. 



Alternatively the factor 3 could be included in the definition of the tensor operator, as sometimes 
done in the literature. In practice what it matters is the product CjOj. 
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Taking i c = we reobtained the corresponding expressions of Ref . [ 03 : • Details will be 
given elsewhere. 

4. Fit and conclusions 

In Table 0] we present the masses of the resonances which we have interpreted as be- 
longing to the [70, + ] or to the [70, 2 + ] multiplet. For simplicity, mixing of multiplets is 
neglected in this first attempt. The resonances shown in column 8, correspond to either 
three stars ("very likely") or to two stars ("fair") or to one star ("poor") status, according 
to Particle Data Group (PDG) [E2]- Therefore we used the full listings to determine a 
mass average in each case. The experimental error to the mass was calculated as the 
quadrature of two uncorrelated errors, one being the average error from the same refer- 
ences and the other was the difference between the average mass and the farthest observed 
mass. For the Pn(2100)* resonance we report results from fitting the experimental value 
of Ref. [EH], as being more recent than the average over the PDG values. Note that the 
observed mass of Ref. [ EH! i s m agreement with the recent coupled channel analysis of 
Manley and Saleski [ I34j . 

Several remarks are in order. Due to its large error in the mass, the resonance Fi5(2000) 
could be either described by the | 2 iV[70, 2+]5/2+) state or by the | 4 iV[70, 2+]5/2+) state 
(inasmuch as they appear separated by about 60-70 MeV only, in quark model studies, 
see, e.g., [ESI ESI)- Here we identified F 15 (2000) with the | 4 iV[70, 2+]5/2+) state because 
it gives a better fit. Regarding the F 35 (1905) resonance there is also some ambiguity. In 
Ref. [EH] it was identified as a | 4 A[56, 2+]5/2+) state following Ref. [ES], but in Ref. 
[122] the interpretation | 4 A[70, 2 + ]5/2 + ) was preferred due to a better x 2 fit and other 
considerations related to the decay width. Here we return to the identification made 
in Ref. [ El] and assign the | 4 A[70, 2 + ]5/2 + ) state to the second resonance from this 
sector, namely -F 35 (2000), as indicated in Table 0] We hope that an analysis based on 
configuration mixing and improved data could better clarify the resonance assignment in 
this sector in the future. Presently the resulting Xdof is about 0.83 and the fitted values 
of q are given in Table [TJ Besides the seven fitted masses Table EJ also contains few 
predictions. 

We found that the contributions of SlS l c and s l S l c are nearly equal when treated as 
independent operators. Therefore, for simplicity, we assumed that they have the same 
coefficient in the mass operator. One can see that the spin-spin interaction given by O4 
is the dominant interaction, as in the [56, 2 + ] multiplet [ 19] or in the [56, 4 + ] multiplet 
[ I22j. Thus the main contributions to the mass come from 0\ and O4. It is remarkable 
that Ci and c 4 of the multiplets [56, 2 + ] and [70, £ + ], both located in the N = 2 band, 
are very close to each other. In terms of the present notation the result of Ref. [ for 
[56, 2+] is ci = 541 ± 4 MeV and c 4 = 241 ± 14 MeV as compared to a = 555 ± 11 MeV 
and C4 = 261 ± 47 MeV here. Such similarity gives confidence in the large N c approach 
and in the present fit. In Fig. 1 the presently known values of ci,C2 and C4 with error 
bars are represented for the excited bands studied within the large N c expansion: N = 1 
is from Ref. [E3> N = 2 from Ref. [QHl and from the present work, N = 4 from Ref. [ 
122]- Extrapolation to higher energies, N > 4, suggests that the contribution of the spin 
dependent operators would vanish, while the linear term in N c , which in a quark model 



Table 4: The partial contribution and the total mass (MeV) predicted by the 1/N C expansion as compared with the empirically 
known masses. 



1/N C expansion results 







Partial contribut 


ion (MeV) 


Total (MeV) 


Empirical (MeV) 


Name, status 






ciOi 


c 2 2 


c 3 3 


C404 








4 iV[70,2+ 


12+ 

1 2 


1665 


31 


42 


217 


1956 ± 95 


2016 ± 104 


Fi 7 (1990)** 


2 iV[70,2+ 


15 + 
1 2 


1665 


10 





43 


1719 ±34 






4 iV[70,2+ 


15 + 
1 2 


1665 


-5 


-106 


217 


1771 ± 88 


1981 ±200 


F 15 (2000)** 


4 iv[7o,o+; 


1 2 


1665 








217 


1883 ± 17 


1879 ± 17 


P 13 (1900)** 


2 iV[70,2+ 


13 + 
1 2 


1665 


-16 





43 


1693 ± 42 






4 iV[70,2+ 


1 2 


1665 


-31 





217 


1851 ±69 






2 iV[70,0+ 


1 2 


1665 








43 


1709 ± 25 


1710 ±30 


Pn(mo)*** 


4 iV[70,2+ 


1 2 


1665 


-47 


149 


217 


1985 ± 26 


1986 ± 26 


Pn(2100)* 


2 A[70,2+j 


5 + 
2 


1665 


-10 





87 


1742 ± 29 


1976 ± 237 


P 35 (2000)** 


2 A[70,2+] 


3 + 
2 


1665 


16 





87 


1768 ± 38 






2 A[70,0+] 


1 + 
2 


1665 








87 


1752 ± 19 


1744 ± 36 


P 3 i(1750)* 



00 
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Figure 1. The coefficients Cj vs N from various sources: for N — 1 from Ref. [E] for 
N = 2 from Ref. [El] (lower values) and present work (upper values), for N = A from 
Ref. [1221- The straight lines are to guide the eye. 
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picture would contain the free mass term, the kinetic and the confinement energy, would 
carry the entire excitation. Such a behaviour gives a deeper insight into the large N c mass 
operator and is consistent with the intuitive picture developed in Ref. [EZj where at high 
energies the spin dependent interactions vanish as a consequence of the chiral symmetry 
restoration. 

Finally note that the contributions of O2 and O3 lead to large errors in the coefficients 
Q obtained in the \ 2 fit, which could possibly be removed with better data. The operator 
O3 containing the tensor term plays an important role in the reduction of xlof an d it 
should be further investigated. 

In conclusion we believe that the mass fit of the [70, £ + ] baryons is encouraging in the 
present form. It is consistent with the conclusions of previous studies that the dominant 
interaction is of a spin-spin nature. It would be interesting to extend this analysis such as 
to incorporate configuration mixing with the [56, 2 + ] states. In this way one could obtain 
a complete description of the N = 2 band baryons, leaving alone the [20, 1 + ] multiplet 
for which no candidate has been found so far. On the other hand improved experimental 
data is highly desirable. 
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Appendix 

Here we express a wave function of a given symmetry described by the partition [/] 
by using fractional parentage coefficients. They are related to the isoscalar factors of the 
permutation group [l2*7j. For one-body operators we need one-body fractional parentage 
coefficients. In this way one can decouple the last particle from the rest. In the simple 
case where the spatial wave function contains only one excited quark, for example having 
the structure (0s) Arc_1 (0c?) (two units of orbital excitation), and symmetry [N c — 1, 1] one 
can show that 



|[JV c -l,l](0 S )^- 1 (0d),2+> = ^^1*^,(0^^- VwCOd) 

-^i*^ ((0sfc- 2 (0d)) ^(Oa) • (14) 

In the case of the spin-orbit operator one can see that only the first term contributes (the 
operator acts on the last particle only). Its matrix element is proportional to the square 
of the coefficient of the first term, i. e. with N ^~ 1 which for large N c gives to the spin-orbit 
the order 0{1) in powers of 1/N C . 
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On the other hand for a symmetric state [N c ] with the same structure one obtains 
\[N c }(0sf^(0d)X) = t /^[;v c -i](0^-V[i](0d) 




which implies that the spin-orbit operator is of order 0(1/ N c ). Both results are in agree- 
ment with Eq. (j2J. 
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